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Abstract 

We study the effect of radiation damping on the classical scattering of 
charged particles. Using a perturbation method based on the Runge-Lenz 
vector, we calculate radiative corrections to the Rutherford cross section, 
and the corresponding energy and angular momentum losses. 

1 Introduction 

The reaction of a classical point charge to its own radiation was first discussed 
by Lorentz and Abraham more than one hundred years ago, and never stopped 
being a source of controversy and fascination Nowadays, it is proba- 

bly fair to say that the most disputable aspects of the Lorentz- Abraham-Dirac 
theory, like runaway solutions and preacceleration, have been adequately un- 
derstood and treated in terms of finite-size effects (for a review see Ref. ;4t). In 
any case, radiation damping considerably complicates the equations of motion 
of charged particles, and for many basic problems, like Rutherford scattering, 
only numerical calculations of the trajectories are available |nilE|- In this paper 
we study the effect of radiation reaction on the classical two-body scattering 
of charged particles. Following Landau and Lifshitz 0, we expand the elec- 
tromagnetic force in powers of (c is the speed of light), up to the order 
where radiation damping appears. Then, using a perturbation technique 
based on the Runge-Lenz vector [Jj , we calculate the radiation damping correc- 
tions to the Rutherford deflection function and scattering cross section, and the 
corresponding expressions for the angular momentum and energy losses. 

This paper is organized as follows. In Sec. 2 we obtain the radiation damping 
force on a system of charged particles, from the expansion of the electromagnetic 
field in powers of 1/c. The equations of motion for a two-body system with 
radiation reaction are discussed in Sec. 3, and in Sec. 4 we use the Runge-Lenz 
vector to calculate the radiation effect on classical Rutherford scattering. Some 
final remarks are made in Sec. 5. 
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2 The radiation damping force 

In this section we reproduce, for completeness, the derivation of the radiation 
damping force given by Landau and Lifshitz 7 . We start from the electromag- 
netic potentials (f>{r,t) and A{r,t), created by the charge and current densities 
p(r, t) and J(r, t), 

^ir,t) ^ I ^d^r' , (1) 

Here, i? = |r — r'| and tj^ = t — R/c is the retarded time. The electric and 
magnetic fields, E and B, are obtained from the potentials as 

E = -V<^(r,i)-i^:^ , B = S7xAir,t). (3) 
c ot 

We want to calculate the electromagnetic force on a charge q, 

F = gE + -V X B, (4) 

c 

as a series in powers of 1/c. In order to do this, we expand p(r', tn) and J(r', tji) 
in Taylor series around = 

p{r',tn) = p{r',t)+ ' ' I ^ ■ i-y ^ ) I 



dt \ c 2 dt^ 



1 d^p{r\t) f R 
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J(r',(«) = J(rM) + ^^^(-^)+0(<:-'). (6) 

Substituting these expansions in Eqs.(^) and ||2J), and noting the charge conser- 
vation relation, 

d 



dt 



J pir\t)dh' = 0, (7) 



we obtain 



lAir,t) . l,J^^d^r'-^§-Jjir',t)d^r' + Oic-^). (9) 
With the gauge transformation 



r,i ^ 0r,t -— 

c ot 
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A{r, t) 



^(r,i)+Vx(r, t) 



(10) 



where 



^^'■''^^dV + 0(c-4) 



-A{v,t) 



1 /■ J(r'>i)^3„/ 



(11) 



(12) 



we can rewrite Eqs. (jHl and © as 

</.(r,t) = 

■ 1 J(r',t)dV-^^ I Rp(r',t)dV + 0(c-4)(13) 
For a set of point charges qu , with positions (t) and velocities (t) , we have 

p{r,t) = ^gfc^(r-rfe(t)), (14) 

fe 

J(r,t) = ^gfeVfe(t)5(r-rfc(t)) , (15) 
fe 

and the potentials become 



9fe 



i?fe(0 



+ 0(c-4) 



(16) 



1 



1 d Rfe(t) 



fe ' ' fe 



fe fe 



with Rfe (<) = r — r^: (t) . Carrying out the time derivatives in Eq. H17|l we obtain 



-A{r,t) 
c 



-y 

2c2 ^ 



gfeVfe(t) gfcRfe(t).Vfe(t) 



i?fe(t) 



Rfe(i) 



3c= 



■E9feafeW + C(c 



(18) 



where a/j(i) is the acceleration of particle k. 

The potential given in Eq. H16|l accounts for the Coulomb interaction. The 
first term in Eq. (|18() . of order l/c^, introduces magnetic and retardation effects, 
and can be used to set up the Darwin lagrangian 2 . The last term in Eq. (jl8|l . 
of order 1/c^, gives the radiation damping electric field 



2 \ - dak 



(19) 
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and a null magnetic field {A is independent of r in this order). Introducing 
the electric dipole of the system, D = J2k1k''^k, the radiation damping field of 
Eq. |{T^ can be written as 

showing that it represents the reaction to the electric dipole radiation emitted 
by the whole system. 

The radiation damping force on charge qi is then 

T^(.i) T-. 2 \ - dak 

K1 = Q^Kd = ^2^ l^lk^ ■ (21) 
fc 

It should be stressed that radiation reaction is not just a self-force — it gets 
contributions from every particle in the system. Only for a single accelerating 
charge q the radiation damping force reduces to the Abraham-Lorcntz self- 
interaction 

2 q2 da 



3 Two-body motion with radiation damping 

Let us consider a system of two charged particles. Taking radiation damping 
into account, their equations of motion read 

d\i qiq2 V 2 qi d 

d^V2 qiq2 Y 2 q2 d 

where t — Yi r2 and is the mass of particle i. In these equations we have 
discarded the c^^ terms that account for the variation of mass with velocity and 
the Darwin magnetic and retardation effects. These terms do not interfere with 
our treatment of radiation damping, and their effect on Rutherford scattering 
is discussed in Refs. [71|S]. 

Subtracting Eq. from (jJU we find 

d^Y qiq2 Y 2 f qi q2 \ d 

-{qiai + q2Si2) , (25) 



dt^ fi 3c^ \ mi m2 / dt 

where fi = mim2 / {mi + 7712) is the reduced mass. From equations H23|) . 124|l 
and (|25|1 . it is easily shown that, keeping only the lowest order (c°) terms, 

f qi q2\ d\ 
qia.1 + (7232 = A* Hrj ■ (26) 



(27) 



^ mi m2 / dt'^ 

Substituting this result in Eq. H25|l we obtain 

d^Y qiq2 Y 2(f' d^Y 
dt"^ fi 3fic^ dP ' 
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where 

^ = ,(^-^). (28) 

\mi 7712/ 

In the fixed target hmit, 7772 — ^ 00, Eq. I|27|l becomes the nonrelativistic 
Lorentz- Abraham equation of motion. It is interesting to see that two-body 
recoil effects appear in Eq. H27|l not only through the reduced mass /i, but also 
via the effective charge q. In particular, if qi/mi = (72/^2 we have q — 0, and 
there is no radiation reaction even though both particles are accelerating. This 
is related to the fact that, in this case, there is no electric dipole radiation from 
the system. 



4 Radiative correction to Rutherford scattering 



In the absence of perturbations, Rutherford scattering conserves the total energy 
E = \lJiv^ + qiq2/r, the angular momentimr L = /ir x v, and the Runge-Lenz 
vector 

A = V X L + qiq2 f . (29) 

Here, v — dr/dt is the relative velocity and f — r/r is the radial unit vector. 
These conserved quantities are not independent: it is easily seen that A • L = 
and 

2ELyfi + {q.q^f = {vaLf + {q.q^f , (30) 

where vq is the initial (asymptotic) velocity. Taking the scalar product r • A, 
one finds the Rutherford scattering orbit 



r{ip) 



A cos 1^ - qiq2 ' 



(31) 



where is the angle between r and A. During the collision, Lp changes from 
—ipQ to (po, where 



ipQ = cos '^{qiq2/A) = tan ^{vQL/qiq2) 



(32) 



The scattering angle is = tt — 2(^0, and from Eq. H32|l we obtain the Rutherford 
deflection function 

e{L) = 2i&n-\qiq2lvoL). (33) 

Note that for charges of the same sign the scattering angle is positive, and for 
opposite charges 6 is negative (we take L and vq as always positive). 

When radiation damping is considered, i?, L and A are no longer conserved. 
In particular, from Eq. (|27() we can show that the Runge-Lenz vector changes 
at the rate 



X L + V X r X 



dA _ 2g2 r 1 d^Y 
'dt 

The total change of A during the collision is then 
2g2 



dt^ 



SA 



3c3 



dt 



1 d^r / d^r 

--rr X L + V X r X — — 

^ dt^ V dt^ 



(34) 



(35) 
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The change of the Runge-Lenz vector is of order c~^. Keeping the same order 
of approximation, we can substitute in the integrand of Eq. (|35|) the results of 
unperturbed Rutherford scattering. We obtain 



<5A 



dt 



A - qiq2 f 



(36) 



which is further simplified by a change of variable from time t to angle if. Still 
working to order c~^, we have 



dt — f^—dtp , 



and 



6A 



3c3 jiL 



dip 



A - qiq2 f 



Substituting r{ip) from Eq. (|^ . the above integral reduces to 



3c3 L^A 



dip {A cos (p - qiq2) (A - qiq2 cos p) 



which is easily calculated. Using Eq. H32|l . the result is written as 



SA 



W glg2«0 

3c3 £2 



1 + (woi/9192) 



j9i92 



tan ^(uoi^/9192) 



(37) 
(38) 

(39) 
(40) 



According to Eq. (|32|) , the change in the Runge-Lenz vector modifies the asymp- 
totic angle pq by 

X '?i92 5 A 



and the scattering angle 9 by (see Ref. [Z]) 

56 = —SpQ . 
The deflection function is then given as 

0{L) = 2tSin-\qiq2/voL) + Se{L) 



(42) 



(43) 



where the first term is the Rutherford relation, and the radiation damping 
correction is 



29^(9192) 
3c3 



1 



1 + {vQL/qiq2)^ vqL 



o9i92 , -1, J- I . 
-3 — ^tan (woi/9192) 



(44) 



From these equations we can also obtain LiQ). To order c the result is 



L(d) = ^cot(0/2) 
Wo 



9^ /wo\3 



9i92 ^ c 



-j 



(45) 
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Figure 1: Angular dependence of the radiative correction to Rutherford's de- 
flection function. Positive (negative) angles correspond to the scattering of like 
(unlike) charges. 



where 



1 sin^{e/2) 



[(5 - cos e) cot(6l/2) - 3(7r - 9)] 



(46) 



6cos5(6i/2) 

A plot of X{9) is shown in Fig. ^ As already mentioned, positive angles are 
reached by like-sign charges, and negative angles by oppositely charged particles. 
We see that the radiative correction is limited if the Coulomb force is repulsive, 
and is strongly divergent for backscattering {9 — tt) in an attractive Coulomb 
field. 

The scattering cross section can be calculated from the deflection as 



da 



1 



L dL 



sine d9 



where p = ^vq is the initial momentum. With Eqs. H45|l and H46|l we get 



da 

dh 



dn 



where 



dn 
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c ' 



1 



2^iVQ^J sin'' (6*72) 
is the nonrelativistic Rutherford cross section, and 



_ 1 s\v?{9/2) 
^(^)-2cos5(^?/2) ^^^-^^^^ 



COS! 



3 sin 61] 



(47) 



(48) 



(49) 



(50) 
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Figure 2: Angular dependence of the radiative correction to the Rutherford 
cross section. Positive/negative angles are the same as in Fig.^ 



The function ^(6*) is shown in Fig. [21 At large angles, close to backscattering, 
£^(9) has the limits 



-96tt{9 + tt) 



(9 ^ -tt) 



(51) 
(52) 



The angular momentum loss (or gain) can be calculated with similar meth- 
ods. With radiation damping, the time derivative of L is given by 



dL 

'dt 



2q^ (Py 



(53) 



which, integrated on the unperturbed Rutherford trajectory, gives the total 
change of angular momentum in the scattering process, 



6L = 



4(?2 qiq2VQ 



1 arctan 



vqL 



3c3 L2 

At a given scattering angle, the angular momentum change is 

4 /llr, \ 3 



3 qiq2 



where 



X{9) = tan2(0/2) 



n -9 



tan(6'/2) 



(54) 

(55) 
(56) 
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Figure 3: Angular dependence of the change in angular momentum. Posi- 
tive/negative angles are the same as in Fig. Q] 



This function is shown in Fig.|21 

The energy loss is readily calculated by differentiating Eq. H30fl . 



5E 



Inserting the expressions for 6K. and (5L we obtain 



or, in terms of the scattering angle 

5E 

"e 



1 



f <M2 

\voL 



arctan 

\q1Q2J 



4 q 



^2 



m 



3 qiq2 V c 

where ^ (0) is the same function given in Eq. H50() and shown in Fig. |2] 



(57) 



(58) 



(59) 



5 Final comments 

Our discussion of radiation damping corrections to Rutherford scattering ig- 
nored relativistic effects like retardation, magnetic forces, and the mass- velocity 
dependence. These effects give contributions of order to the deflection func- 
tion and cross section (see Ref. 0), and are generally more important than 



9 
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Figure 4: Rutherford scattering to order c^^. The projectile velocity is 0.4 c, and 
the target has infinite mass. The two electric charges are of the same magnitude, 
like (unlike) signs corresponding to positive (negative) scattering angles. The 
nonrelativistic Rutherford cross section is given by the dotted lines. The dashed 
lines incorporate c"^ corrections, and the solid lines include the c^^ radiation 
damping effects. 



the c^^ radiative corrections we have obtained. They were not considered here 
because, as already mentioned, this would not change our results: a cor- 
rection to the nonrelativistic Rutherford trajectory only adds terms to our 
perturbative calculation of radiation damping. We can easily write the complete 
(up to c~^) expansion of the deflection function and scattering cross section by 
putting together the results of Ref. and the present paper. For example, the 



differential cross section to order 



reads 



da 
where 



c ' 
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h{9) = - tan2(6l/2) [1 + (tt - 9) cot 9] + 1 



(60) 



(61) 



and AI = mi + m2. As discussed in 7 , the first corrective term accounts 
for the variation of mass with velocity, and the second includes magnetic and 
retardation effects. The last one is the radiative correction calculated in the 
previous section. It is interesting to note that magnetic and retardation effects 
simply renormalize the cross section by an angle independent factor. 

In Fig. 0]we show the differential cross section for the scattering of a charged 
particle with vq — 0.4 c on a fixed target, of equal {9 > 0) or opposite (9 < 0) 
charge. The dotted lines give the nonrelativistic cross section, and the dashed 
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ones show the effect of the relativistic mass correction (retardation and 
magnetic forces do not show up on a fixed target). The sohd fines bring in the 
radiation damping effect, as given in Eq. (|(jO|l . We see in Fig. 21 that radiation 
damping has a very small effect when the charges repel each other. But for 
an attractive Coulomb force the radiative correction is quite important (as also 
seen in Fig. |2Jl, creating a plateau- like structure in the angular distribution. 
Even though our perturbative results are not reliable for large corrections, such 
structure is very similar to what is found in "exact" numerical calculations 

A final point we wish to comment on is why our results are not plagued 
by runaway solutions. The reason is that the Runge-Lenz based perturbative 
calculation presented here follows essentially a "reduction of order" approach, 
such as described in Refs. 12 ^]. This effectively eliminates the additional 
degrees of freedom introduced in the equations of motion by the time derivative 
of acceleration, yielding only physically acceptable solutions. 
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